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Abstract 



We present a novel scheme for nuclear structure calculations based on realistic nucleon- 
nucleon potentials. The essential ingredient is the explicit treatment of the dominant inter- 
action-induced correlations by means of the Unitary Correlation Operator Method (UCOM). 
Short-range central and tensor correlations are imprinted into simple, uncorrected many- 
body states through a state-independent unitary transformation. Applying the unitary trans- 
formation to the realistic Hamiltonian leads to a correlated, low-momentum interaction, 
well suited for all kinds of many -body models, e.g., Hartree-Fock or shell-model. We em- 
ploy the correlated interaction, supplemented by a phenomenological correction to account 
for genuine three-body forces, in the framework of variational calculations with antisym- 
metrised Gaussian trial states (Fermionic Molecular Dynamics). Ground state properties 
of nuclei up to mass numbers A < 60 are discussed. Binding energies, charge radii, and 
charge distributions are in good agreement with experimental data. We perform angular 
momentum projections of the intrinsically deformed variational states to extract rotational 
spectra. 
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1 Introduction 

The advent of realistic nucleon-nucleon (NN) potentials has created a supreme 
challenge and opportunity for nuclear structure theory: the ab initio description 
of nuclei. Several families of realistic or modern NN potentials have been devel- 
oped over the past three decades — among others the Argonne, Bonn, and Ni- 
jmegen potentials [1,2,3,4]. The latest members of these families reproduce the 
NN-scattering data and the deuteron properties with high precision. All realistic 
potentials exhibit a quite complicated operator structure with substantial contribu- 
tions from non-central terms, most notably the tensor and spin-orbit interaction. 
Besides these, explicit momentum dependent terms or quadratic orbital angular 
momentum and spin-orbit operators enter. The most recent potentials also employ 
charge- symmetry and charge-independence breaking terms to further improve the 
agreement with the experimental phase- shifts. 

Given a realistic NN-potential we are faced with the difficult task of solving the 
nuclear many-body problem. Ideally, one would like to treat the (non-relativistic) 
quantum many-body problem ab initio, i.e., without introducing further approx- 
imations. So far, ab initio solutions utilising realistic potentials are restricted to 
light nuclei with A < 12. Extensive studies on ground state properties and excita- 
tion spectra in this mass range have been performed using, e.g., Green's Function 
Monte Carlo methods [5,6] and the no-core shell model [7,8]. 

The results of these ab initio calculations clearly show that in addition to the 
realistic NN-potential a three-nucleon force is inevitable to reproduce the experi- 
mental data on light nuclei. The Argonne group has constructed a series of phe- 
nomenological three-nucleon potentials [9] to fit the binding energies and spectra 
of light nuclei to experiment. In this way the results for ground states and low- 
lying excitations are in good agreement with the experimental findings in the whole 
accessible mass range. Recent developments in chiral effective field theories [10] 
provide a framework for a consistent derivation of two-, three- and multi-nucleon 
forces from more fundamental grounds. 

Our aim is to describe the structure of larger nuclei on the basis of realis- 
tic nucleon-nucleon interactions while staying as close as possible to an ab initio 
treatment of the many-body problem. The step towards larger particle numbers re- 
quires a truncation of the full many-body Hilbert space to a simplified subspace 
of tractable size. In the simplest approach, the many-body state is described by a 
single Slater determinant, as, e.g., in the Hartree-Fock approximation. More elab- 
orate approximations, the multi-configuration shell-model for example, allow for 
many-body states which are represented as a superposition of several Slater de- 
terminants. When combining realistic nucleon-nucleon interactions with simple 
many-body states composed of single or few Slater determinants a fundamental 
problem arises. Those states do not allow for an adequate description of the strong 
short-range correlations induced by the realistic NN-potential [1 1,12]. 

Already in the deuteron, structure and origin of these correlations are appar- 
ent. We consider the spin-projected two-body density matrix p" s '_ x M X r ) as function 
of the relative coordinate r = jci - x 2 of the two particles, resulting from an ex- 
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Fig. 1. Spin-projected two-body density pj jL (r) of the deuteron calculated with the AVI 8 
potential. Shown is the iso-density surface for 0.005 fm . 
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Fig. 2. (a) Central part of the AV18 potential in the (5, T) = (0, 1) channel, (b) Two-body 
density distribution p^ 07 . =1 (r) for 4 He resulting from an ansatz wave function using 

a Slater determinant of Gaussian single-particle states ( ) compared to a realistic 

two-body density including the interaction-induced central correlations ( ). 

act calculation with the Argonne V18 potential (AV18) [1,13,12]. Figure 1 depicts 
iso-density cuts of the two-body density for M s = and M s = ±1, respectively, 
which corresponds to parallel and antiparallel alignment of the two nucleon spins. 
Two dominant structural features appear: (?) At small particle distances \r\ the two- 
body density is fully suppressed, i.e., the probability of finding two nucleons closer 

(1) 

than r ~ 0.5 fm is practically zero, (ii) The angular structure of p s=x M (r) de- 
pends strongly on the spin orientation. For parallel spins the probability density is 
concentrated along the quantisation axis ("dumbbell"), for antiparallel spins it is 
constricted around the plane perpendicular to the spin direction ("doughnut"). 

These structures are the manifestation of two types of interaction-induced cor- 
relations which govern the nuclear many-body problem: (i) short-range central cor- 
relations and (ii) tensor correlations. 

The central correlations are induced by the strong repulsive core of the central 
part of the interaction. Figure 2 shows the radial dependence of the central part 
of the AVI 8 potential in a specific spin-isospin channel. Due to the short-range 
repulsion the two-body density, shown for 4 He, is completely suppressed within 
the region of the repulsive core. A typical ansatz for the many-body wave function 
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Fig. 3. Dependence of the tensor interaction energy on the spatial orientation of the two 
spins for fixed antiparallel (a) and parallel spins (b), respectively. The numbers indicate the 
values of Vt for each of the configurations. 



using a Slater determinant of Gaussian (or any other) single-particle states is not 
capable of describing this correlation hole. The dashed line in Fig. 2(b) shows that 
the Slater determinant ansatz leads to a significant probability of finding the nu- 
cleons within the core. Even with a superposition of a moderate number of Slater 
determinants one is not able to describe the strong short-range central correlations 
caused by realistic potentials. 

The origin of the tensor correlations is the long-range tensor interaction, partic- 
ularly in the (S, T) = (1,0) channel. As seen in Fig. 1, the tensor correlations are 
revealed through the dependence of the two-body density on the spatial orientation 
of the two nucleons with respect to their spins. This can be understood from the 
structure of the tensor operator Si 2 = 4( cr i ■ r )( cr 2 ■ r ) _ (p~\ • o"i)- For illustration 
we consider the interaction energy V T = -[4(o"i • r)(cr 2 • r) - (o^ • er 2 )], where the 
minus sign was introduced to account for the sign of the radial dependence of the 
tensor interaction. Spin vectors and relative coordinate enter in the same way as for 
the magnetic dipole-dipole interaction. Figure 3 gives the values of V T for differ- 
ent spatial orientations of two fixed antiparallel and parallel spins, respectively. For 
antiparallel spins, those configuration with r perpendicular to the spin direction are 
energetically favoured. Contrariwise, for parallel spins a fully aligned arrangement 
of spins and relative coordinate is preferred. This explains the structure of the two- 
body density shown in Fig. 1 . It is clear that these tensor correlations between the 
two spins and the relative coordinate cannot be described adequately by a single 
Slater determinant. 

The aim of this paper is to devise and apply a method to include central and 
tensor correlations explicitly into simple model spaces. We discuss this approach, 
the Unitary Correlation Operator Method (UCOM), in Sec. 2. The structure of the 
correlated realistic interaction is presented in Sec. 3 with application to the AVI 8 
potential. Section 4 introduces a simple yet powerful variational model for the 
treatment of the nuclear many-body problem based on the Slater determinant of 
Gaussian wave packets used in Fermionic Molecular Dynamics (FMD). Ground 



state calculations up to mass numbers A < 60 are performed, and binding ener- 
gies, charge radii, and charge distributions are compared to experiment. Finally, 
Sec. 5 discusses the projection onto angular momentum eigenstates of the intrinsi- 
cally deformed variational states and compares the resulting rotational spectra with 
experiment. 

2 Unitary Correlation Operator Method (UCOM) 

2. 1 Concept 

The basic idea of the Unitary Correlation Operator Method (UCOM) is the 
following: Introduce the dominant short-range central and tensor correlations into 
a simple many-body state by means of a state-independent unitary transformation. 
The unitary operator C describing this transformation applied to an uncorrelated 

«<V 1 

many-body state \¥) leads to a correlated state |*F) , 

m = c m . CD 

In the simplest case, the uncorrelated state is a Slater determinant. The correlated 
state |¥), however, is no longer a Slater determinant. Due to the complex short- 
range correlations introduced by the correlation operator C, any expansion of \¥) 
in a basis of Slater determinants will require a huge number of basis states. 

The representation of short-range correlations in terms of a unitary state-inde- 
pendent operator is technically very advantageous. Calculating a matrix element 
(^l O |4") of an operator O with correlated states is equivalent to calculating the 
matrix element (¥1 O |*F') of the correlated operator 

6 = C t OC (2) 

using uncorrelated states. For the treatment of the many-body problem it is gener- 
ally more convenient to correlate all operators of interest and to evaluate expecta- 
tion values or matrix elements with the uncorrelated states. 

In accord with the two types of correlations discussed above, we decompose 
the correlation operator into separate unitary operators C n and C r for tensor and 
central correlations, respectively, 

C = C n C r = exp[ - i 2] gCKij] exp[ - i ^ g rJj \ . (3) 

Each of the correlation operators can be written as an exponential involving a Her- 
mitian generator. Since the correlations considered here are induced by a two-body 
potential, the generators are also assumed to be two-body operators. The detailed 
form of the two-body generators g,. and g n reflects the structure of the central and 



1 Throughout the paper, hats identify correlated quantities; upright symbols indicate oper- 
ators. 



tensor correlations, as discussed in the Sec. 1 . We will construct these generators in 
the following sections. 

2.2 Central Correlations 

The short-range repulsion in the central part of the NN-interaction prevents the 
nucleons in a many-body system from approaching each other closer than the extent 
of the repulsive core, i.e., the two-body density matrix exhibits a correlation hole 
at small interparticle distances. One way to incorporate these correlations into a 
simple many-body state is by shifting each pair of nucleons apart from each other. 
The shift has to be distance-dependent since it should only affect nucleon pairs 
which are closer than the core radius. 

Using this picture, we can construct the following ansatz for the Hermitian gen- 
erator of the corresponding unitary transformation (3): radial shifts are generated 
by the component of the relative momentum q = i[pj - p 2 ] along the distance 
vector r = Xi - x 2 of two particles: 

q, = ^[q-7 + ?-q]- (4) 

To describe the dependence of the shift on the particle distance we introduce a 
function s{r) and define the following Hermitian generator for the central correla- 
tion operator C r [11]: 

g r = -[s(r)q r + q r s(r)] . (5) 

To illustrate the effect of the correlation operator c r = exp[-ig r ] in two-body 
space 2 , we apply it to a two-body state |i/^) = |O cm ) <g> \<p). The correlation oper- 
ator does not act on the centre of mass component |O cm ) by construction. For the 
correlated relative part |0) we obtain in coordinate representation [11]: 

R-(r) 



(r\fo = (r\ c r |0) = — — ^JR r {rj (R-(r^\cf>) . (6) 

r r 

Hence, the application of the correlation operator corresponds to a norm conserving 
coordinate transformation r i-> R-{r) L with respect to the relative coordinate. The 
function R-(r) and its inverse R+(r) are connected to the shift function s(r) by the 
integral equation 



r-R±(r) 



±1, R ± [R*(r)]=r. (7) 



For slowly varying shift functions s(r), the correlation functions are approximately 
given by R±(r) « r + s(r). The R±(r) are determined by an energy minimisation in 
the two-body system, as discussed in Sec. 3.5. 



2 c stands for the correlation operator in two-body space, whereas C indicates the general 
correlation operator in many-body space. In general, small symbols are used for &-body 
operators in fc-body space, capital symbols for operators in many-body space. 



Since the NN-interaction depends strongly on the spin S and isospin T of the 
interacting nucleon pair, we decompose the two-body generators g r into a sum of 
different generators gf T for each of the four different (S, T) channels 

g r = ^gfn sr , (8) 

S,T 

where H ST is a projection operator onto the spin S and isospin T subspace. The 
correlation operator in two-body space decomposes into a sum of independent cor- 
relation operator for the different channels 

c r = exp[ - i J] d T n ST ] = J] ex PHgf] U st = J] c r Tu ST ■ (9) 

S,T S,T S,T 

This form is very convenient, since it allows us to treat the (S, T) channels sepa- 
rately. For the sake of a concise formulation, we will not write out this spin-isospin 
dependence in the following, but add it when needed. 

2.3 Tensor Correlations 

The strong tensor part of the NN-interaction induces subtle correlations be- 
tween the spin of a nucleon pair and their relative spatial orientation. In order to 
describe these correlations, we need to generate a spatial shift perpendicular to the 
radial direction. This is done by the "orbital momentum" operator 

qn = q-?qr = 2^[lxr-rxl]. (10) 

where 1 is the relative orbital angular momentum operator. Radial momentum -q, 
and the orbital momentum q n constitute a special decomposition of the relative 
momentum operator q and generate shifts orthogonal to each other. The complex 
dependence of the shift on the spin orientation is encapsulated in the following 
ansatz for the generator g n [12]: 

3 
gn = -^(r)[(o- r q n )(o- 2 - r) + (ov r)(tr 2 - q n )] = &(r) & 12 (r, q n ) , (11) 

where s 12 (r, q n ) = \\ip-i ■ q n )(o" 2 • r) + (en • r)(cr 2 • q n )]. The two spin operators 
and the relative coordinate r enter in a similar manner like in the tensor operator 
Si 2 , however, one of the coordinate operators is replaced by the orbital momentum 
q n , which generates the transverse shift. The size and the distance-dependence of 
the transverse shift is given by the function #(r) — the counterpart to s(r) for the 
central correlator. The isospin dependence of the tensor correlator is implemented 
in analogy to the spin-isospin dependence of the central correlator (9) through pro- 
jection operators. 

Again we illustrate the effect of the correlation operator Cn = exp[-ig n ] by ap- 
plying it to a two-body wave function in coordinate representation. For simplicity 




Fig. 4. Construction of the deuteron wave function for the AVI 8 potential starting from 
an uncorrected wave function shown in panel (a). Application of the central correlator 
with a correlation function R+(r) shown in (d) leads to the central correlated wave function 
(b). Application of the tensor correlator with &(r) shown in panel (e) produces a D-wave 
admixture and leads to a realistic deuteron wave function depicted in (c). 

we use LS -coupled two-body states \tp; (LS)J). Starting from a pure L = un- 
corrected wave function with S = 1 , T = 0, and J = 1 , for example, the tensor 
correlator generates a superposition of an L = and an L = 2 state 



(r\ c n 10; (01)1) = cos(3 V2 0(r)) <r|0; (01)1) 

+ sin(3 V2 #(r)) <r|0;(21)l> . 



(12) 



The tensor correlation function &(r) determines the amplitude and the radial depen- 
dence of the D-wave admixture. 

To further illustrate the effect of the central and the tensor correlation operator 
on a two-body state, Fig. 4 details the steps from a simplistic ansatz to the exact 
deuteron wave function for the AVI 8 potential. We start with an uncorrelated L = 
wave function, shown in panel (a), which correctly describes the long-range be- 



haviour but does not contain a correlation hole or a D-wave admixture. Applying, 
in a first step, the central correlation operator with the correlation function R + (r) 
depicted in panel (d) leads to the correlated wave function shown in Fig. 4(b). The 
amplitude is shifted out of the core region towards particle distances where the po- 
tential is attractive. To this end, the shift function s(r) « R+(r) - r has to be large 
within the core radius, but has to decrease rapidly outside the core region. In a 
second step, we apply the tensor correlator for a correlation function &(r) given in 
Fig. 4(e). As seen from Eq. (12) this generates a D-wave admixture and leads to 
a fully correlated wave function depicted in panel (c), which is in nice agreement 
with the exact deuteron solution. Note that the shape of the D-wave component is 
determined by the correlation function #(r) — a long-range L = 2 wave function 
requires a tensor correlation function &(r) of long range. 

2.4 Correlated Operators and Cluster Expansion 

The explicit formulation of correlated wave functions for the many -body prob- 
lem becomes technically increasingly complicated, and the equivalent notion of 
correlated operators proves more convenient. The similarity transformation (2) of 
an operator O leads to a correlated operator which contains irreducible contribu- 
tions to all particle numbers. We can formulate a cluster expansion of the correlated 
operator 

6 = c^oc = 6 [1] + 6 [2] + 6 [3] + • • • , (13) 

where O w denotes the irreducible n-body part [11]. When starting with a £-body 
operator, all irreducible contributions [ " ] with n < k vanish. Hence, the unitary 
transformation of a two-body operator — the NN-interaction for example — yields 
a correlated operator containing a two-body contribution, a three-body term, etc. 

The significance of the higher order terms depends on the range of the cen- 
tral and tensor correlations [11,12,14]. If the range of the correlation functions is 
small compared to the mean interparticle distance, then three-body and higher- 
order terms of the cluster expansion are negligible. Discarding these higher-order 
contributions leads to the two-body approximation 

6 C2 = 6 [1] + 6 [2] . (14) 

Already the example discussed in connection with Fig. 4 shows that for the 
central correlations the two-body approximation is well justified. The range of the 
shift function s(r) « R + (r) - r is smaller than the mean interparticle distance at 
saturation density which is about 1.8 fm. One can show by an explicit evaluation 
[14] that higher order contributions due to central correlations are indeed small for 
the nuclear many-body problem. 

The situation is different for tensor correlations. The tensor correlation function 
#(r) needed to generate the D-wave component of the deuteron wave function is 
necessarily long-ranged (see Fig. 4(e)). In a many-nucleon system, however, the 
tensor correlations between two nucleons will not be established up to the same 
large distances as in the deuteron. The other nucleons interfere and inhibit the for- 



mation of long-range tensor correlations and thus lead to an effective screening 
of the long-range tensor correlations. If we were to use the long-range tensor cor- 
relator suitable for the deuteron also in the many-body problem, this screening 
effect would emerge through substantial higher-order contributions of the cluster 
expansion. Since an explicit calculation of the higher-order tensor contributions to 
the cluster expansion is presently not feasible, we will anticipate the screening of 
the long-range tensor correlations by restricting the range of the tensor correlation 
function [12]. In this way, we can improve upon the quality of the computationally 
simple two-body approximation. Possible residual long-range tensor correlations 
that are not represented in the short-range tensor correlator have to be described by 
the many -body states of the model space. We will come back to the construction of 
the optimal correlation functions for a given NN-potential in Sec. 3.5. 



3 Correlated Realistic Interactions 

3.1 Nuclear Hamiltonian 

We are now applying the formalism of the Unitary Correlation Operator Method 
discussed in Sec. 2 to construct the correlated nuclear Hamiltonian in two-body 
approximation. The starting point is an uncorrelated Hamiltonian for the A-body 
system 

H = T+v =E2^ p ? + i>- (i5) 

1=1 " 1>J=1 

consisting of the kinetic energy operator T and a two-body potential V. For the 
latter, we employ realistic NN-potentials from the Bonn or Argonne family of in- 
teractions. Those are given in a closed operator representation facilitating the use 
within the UCOM framework. The Argonne V14 [15] and the charge independent 
terms of the Argonne VI 8 interaction [1] have the following operator structure 



V Argonne = Y ; K r 0") + vf T (l) 1 2 ] U ST 



\rgonne / 



^ (16) 

+ 2j ^ (r) Sl2 + VT® (1 • S) + ^W Q • S ) 2 ] n ^ ' 
T 

where Ust denotes the projection operator onto spin S and isospin T. In the fol- 
lowing, operators n 5 with a single index always refer to a projection in spin-space 
only. The quadratic spin-orbit term can be rewritten 

(l-s) 2 = |l 2 n 1 -i(l-s) + is 12 (l,l), (17) 

where 

s 12 (l, 1) = 3(o-! • l)(cr 2 • 1) - (<r, • <r 2 ) 1 2 . (18) 
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The non-relativistic configuration space versions of the Bonn potentials [2] are 
parametrised using a different set of operators 



VBonn = J] [ V ^« + \{^f T (r) + vf T (r) q 2 }] U ST 

S,T 

+ £ [v' T (r) s 12 + v£(r) (1 • s)] U 1T . 



T 

Instead of the l 2 and (1 • s) 2 terms in (16), the Bonn potentials employ a non-local 
momentum-dependent term involving q 2 . For the following considerations it is con- 
venient to rephrase this contribution in terms of the radial momentum q r and the 
angular momentum 1 

i{qV 2 (r) + v* 2 (r) q 2 } = ±{q 2 v* 2 (r) + v q2 {x) q 2 } + ^-^- 1 2 . (20) 

3.2 Tensor Correlated Hamiltonian 

The first step to construct the correlated Hamiltonian in two-body approxima- 
tion is the application of the tensor correlation operator c n = exp(-ig^) in two- 
body space. A general way to evaluate the similarity transformation is by utilising 
the Baker-Campbell-Hausdorff expansion 

cJjO c n = exp(ig n )0 exp(-ig n ) 

i 2 
= + i[g n ,0] + -[gn,[gn, ()]] + •■• 



^-L^O = exp(L n )0 



n=0 

In the last line we have introduced a compact notation of the iterated commutators 
in terms of powers of the super-operator L n o = i[gn, °] with the generator g n given 
by Eq. (11). Summing up the full expansion formally leads to the exponential of 
the super-operator. 

First we study the transformation of the various terms of the realistic NN- 
potentials (16) and (19). A minimal set of operators we have to consider in order to 
formulate the tensor correlated interaction is {r, q 2 , l 2 , (1 • s), s 12 , Si 2 (l, 1)}. The dis- 
tance operator r commutes with the generator g n , i.e., it is invariant under tensor 
correlations 

c^rc n = r. (22) 

For the radial momentum q 2 , the Baker-Campbell-Hausdorff expansion terminates 
after the second order and we obtain a closed expression for the tensor correlated 
operator [12], 

c^q 2 c n = q 2 - [#'(r) q r + q,#'(r)] s 12 (r, q n ) + [#'(r) s 12 (r, q n )] 2 , (23) 

which can be further simplified using the identity s J2 (r,q n ) 2 = 9[s 2 + 3(l-s) + (l-s) 2 ]. 
All other basic operators require the evaluation of the full commutator expansion 
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(21). At first order, the following commutators appear [16]: 

[go, si 2 ] = i0(r)[-24IIi - 18 (1 • s) + 3 s 12 ] 
[go. • s)] = i#(r)[-s 12 (q n , q n )] 
[gn,l 2 ]=i^(r)[2s 12 (q n ,q n )] 
[gn,s 12 (l,l)] = i#(r)[7 s 12 (q n , q n )] , 

where we have used the short-hand notation 

si2(qn, qn) = 2r 2 s 12 (q n , q n ) + s i2 (l, 1) - \ s 12 . (25) 

In addition to the original set of operators, a new tensor containing two orbital 
momentum operators is generated 

s i2(q n > qn) = 3(0"! • q n )(o"2 • q n ) - (°"i • °"2) q^ • (26) 

For the calculation of the second order of the expansion (21), the commutator of gn 

and Si 2 (q n , q n ) is required: 

[gn,s 12 (q n ,qo)] = i#(r)[(108 + 961 ^ + (153 + 361 2 )(1 • s) + 15s 12 (l,l)] . (27) 

Again a new operator, 1 2 (l-s), emerges, whose commutator with the generator enters 
into the third order of the Baker-Campbell-Hausdorff expansion: 

[g n ,l 2 (l • s)] = i#(r)[-3 s 12 (q n ,q n ) - {l 2 s 12 (q n ,q n )} H ] . (28) 

This in turn leads to the new operator {1 2 s 12 (q n , q n )} H = 5(1 2 Si 2 (q n , q n )+s 12 (q n , q n )l 2 ), 
whose commutator appears in the fourth order of the expansion: 



[gn,{l 2 Si 2 (q n ,qn)} H ] =i^r)[324n, + 477(1 • s) + 6001 2 + 51s 12 (l,l) 

+ 4771 2 (1 • s) + 1441 4 + 271 2 s 12 (l,l) + 361 4 (1 • s)] . 



Evidently, with increasing order in the Baker-Campbell-Hausdorff expansion, op- 
erators containing successively higher powers of the angular momentum 1 are gen- 
erated. The last three terms of (29) are of fourth order in 1 already. In the following, 
we will neglect the contributions beyond the third order in angular momentum. 
Thus we achieve a closure of the operator set contributing to the Baker-Campbell- 
Hausdorff expansion: 

{r,q 2 ,l 2 ,(l-s),si 2 ,si 2 (l,l), 

Si2(q n > q n X q r si 2 (r, q n ), {1 2 Si 2 (q n , q n )} H , 1 2 (1 • s)} . 

We can represent the super-operator L n defined in (21) as a matrix acting on the 
vector (30) of operators. The summation of the full Baker-Campbell-Hausdorff ex- 
pansion is then reduced to calculating a matrix exponential for the super-operator. 
In this way a closed operator representation of the tensor correlated potential is 
constructed. 
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Fig. 5. Radial dependency of the tensor part of the AVI 8 potential ( ) and the radial 

dependencies ( ) of the tensor (a), central (b), and spin-orbit contribution (c) that 

result from applying the tensor correlator to the isolated tensor part of the interaction. All 
potentials are for the (S,T) = (1,0) channel, and the appropriate tensor correlation function 
is given in Sec. 3.5 (a-correlators). 

The set (30) already contains the operators relevant for the transformation of 
the kinetic energy. We can decompose the kinetic energy T in two-body space into 
a relative, t re i, and a centre of mass contribution, t cm . The latter is not affected by 
the correlations. The relative contribution is further decomposed into a radial and 
an angular term 

1 l 2 

T = ton + trel = U + (q? + -r) . (31) 

Thus the transformation of the kinetic energy can be inferred from the transforma- 
tion properties of the first three operators in the set (30). 

The effect of the tensor correlator on the tensor part of the interaction is illus- 
trated in Fig. 5. For this example we take the isolated tensor part of the AVI 8 poten- 
tial in the (S, T) = (1,0) channel and apply the tensor correlation operator. Through 
the correlation procedure, i.e., the application of the matrix exponential according 
to (21), contributions to the other operator channels (30) are generated. Most no- 
tably, additional central and spin-orbit contributions emerge, which are depicted in 
panels (b) and (c) of Fig. 5. Both, the induced central and spin-orbit contributions 
are attractive at intermediate ranges. Hence, through the correlation procedure, part 
of the tensor attraction is transfered to other operator channels, where it can be 
exploited using simple uncorrelated many-body states. The marginal repulsion at 
short ranges will be completely removed by the central correlation operator as dis- 
cussed in the following. 

3.3 Central and Tensor Correlated Hamiltonian 

The subsequent application of the central correlator to the tensor correlated 
terms of Hamiltonian is technically simpler. In Section 2.2 we have shown that the 
central correlator acts like a norm conserving coordinate transformation when ap- 
plied to a relative two-body wave function (6). Using the transformation properties 
of the wave function in coordinate representation we can immediately derive ex- 
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pressions for the central-correlated operators in two-body approximation [11]. The 
evaluation of the Baker-Campbell-Hausdorff expansion for the central correlations 
is therefore not required. 

The application of the central correlation operator c,. to the distance operator r 
confirms the picture of a coordinate transformation 

clrc r = R + (r), (32) 

where R+(r) is the correlation function introduced in Sec. 2.2. Wherever the dis- 
tance operator r appears in the Hamiltonian, it has to be replaced by a transformed 
distance operator R+(r). This affects, most notably, all radial dependencies of the 
different contributions of the interaction 

cjv(r)c r = v(i? + (r)). (33) 

The components of the relative momentum operator q are also influenced by the 
central correlations. Using the correlated two-body wave function (6), one can show 

that 

t 1 1 t r 

c;q r c r = q r , c;q n c r = — -q n . (34) 

V^(rj ^JRW) *+< r ) 

For the quadratic radial momentum q 2 , which appears in the tensor correlated po- 
tential as well as in the kinetic energy (31), one obtains 

,t 2 . l l l 2 . _ 2 1 1 . 7[i?;'(r)] 2 K'{r) 



CArCr 2 1 [i?' + (r)]2 q '" + q; [R> + (r)f J + 4|&(r)]« 2[R' + (r)¥ ' (35) 

Notice that the transformation of q 2 generates a local potential in addition to the 
momentum-dependent term. 

All basic operators that act only on the angular part of the two-body wave func- 
tion, are invariant under central correlations, e.g., 

c!lc r = l, c!-c r = - . (36) 

r r 

Utilising the unitarity of the correlation operator, one can easily deduce from these 
basic identities the correlated expressions for the composite tensor and spin orbit 
operators contained in the operator set (30). One finds that l 2 , (1 • s), Si 2 , s 12 (l, 1), 
Si2(Qn' Qn)' (1 2 Si2(Qn> Qn)}^ and 1 2 (1 • s) are invariant under similarity transforma- 
tion with the central correlation operator. 

3.4 Properties of the Correlated Interaction Vucom 

Combining the different central and tensor correlated operators discussed in 
the previous sections, we can formulate the correlated many-body Hamiltonian in 
two-body approximation 

H C2 = f [1] + f [2] + V [2] = T + Vucom • (37) 
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(38) 



Its one-body contribution is just the uncorrected kinetic energy T [1] = T. Two- 
body contributions arise from the correlated kinetic energy T [2] and the correlated 
potential V [2] . Together they constitute the correlated interaction Vucom, which is 
the basis for the following nuclear structure studies. 

The generic operator structure of the correlated interaction V UCO m is more com- 
plicated than that of the bare potentials 

vucom = Yj KrW + &7fc)£h + vfr(r)l 2 ] Ust 

S,T 

+ J] [vr(r) (1 • s) + v' T (r) s 12 + if (r) s 12 (l, 1) 

T 

+ v? q (r) s 12 (q n , q n ) + vf qq (r) {l 2 Si 2 (q n , q n )} H 
+ {vf tr \r) q 2 s 12 (r, q n )} H + vf (r) 1 2 (1 • s)] n ir , 

where we use the short-hand notation {• • -} H to indicate explicit Hermitisation. The 
new radial dependencies v ST (r) result from the correlation procedure, i.e., from the 
matrix exponential for the tensor correlator and the coordinate transformation for 
the central correlator. They depend on the radial dependencies of the bare potential 
and on the correlation functions. The structure of the different radial dependencies 
was discussed in Refs. [11,12]. In summary, the two prime effects of the unitary 
transformation are: (i) The short-range repulsive core of the central interaction is 
removed by the central correlator and an effective repulsion in the momentum- 
dependent terms is generated, (ii) Additional attractive central and spin-orbit con- 
tributions as well as new tensorial terms are created by the tensor correlator out of 
the bare tensor interaction. 

Before entering into concrete many-body calculations, we summarise a few key 
properties of the correlated interaction Vucom- First, the correlated interaction is 
phase-shift equivalent to the uncorrected potential by construction. This is a direct 
consequence of the finite range of the correlation functions. The asymptotics of a 
scattering wave function is not altered by the correlation operators and the phase 
shifts are preserved. Hence, the unitary correlation operator provides a tool to gen- 
erate an infinite manifold of phase- shift equivalent NN-potentials originating from 
a single realistic interaction. In addition to the two-body potential the correlation 
operator generates a three-body (and higher-order) interaction which, of course, 
depends on the particular correlation functions used. 

Furthermore, there is an interesting connection to the low-momentum interac- 
tion Viow/t determined by means of renormalisation group techniques [17]. The 
momentum space matrix elements of Vucom are m very good agreement with the 
V\ ^-u matrix elements [12]. Although both approaches are formally quite differ- 
ent, the underlying physics is the same: the high-momentum components of the 
interaction are treated explicitly — by the unitary correlation operator or through 
the renormalisation group procedure — leaving an effective interaction adapted 
to low-momentum model spaces. One major practical advantage of the correlated 
interaction Vucom is that it is given in a closed operator representation (38). De- 
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pending on the particular application, one can easily compute the relevant matrix 
elements, e.g., in a plane wave, oscillator, or non-orthogonal Gaussian basis. The 
renormalisation group method just provides numerical values for the momentum- 
space matrix elements. 



3.5 Optimal Correlation Functions 

A crucial step is the construction of optimal correlation functions for the use 
in many-body calculations. These correlation functions depend, of course, on the 
bare potential, but they should not depend on the nucleus under consideration. The 
aim is to fix a set of optimal state-independent correlation functions which define 
a fixed correlated interaction Vucom- When constructing the correlation functions 
one, therefore, has to disentangle state-dependent properties and state-independent 
features. 

This is an issue in particular for the tensor correlations in the (S, T) = (1,0) 
channel. The long-range tensor correlation function &(r), which was used in Fig. 
4(e) to reproduce the D-wave admixture of the exact deuteron solution, is certainly 
not appropriate for a many-body system. As already discussed in Sec. 2.4, long- 
range tensor correlations of a pair of nucleons are screened through tensor interac- 
tions with other nucleons. One way to isolate the essential state-independent con- 
tributions is to restrict the tensor correlation functions to short ranges. Long-range 
tensor correlations, which are not accounted for explicitly by the short-range corre- 
lator, have to be described through the degrees of freedom of the many-body states. 
For the central correlations this restriction to short ranges emerges automatically, 
since the strong repulsive core and the induced correlation hole are short-ranged 
themselves. 

Another motivation to consider only correlation functions of short range is the 
validity of the two-body approximation. The two-body approximation is appropri- 
ate as long as the correlation range is small compared to the mean interparticle 
distance. In this case, the probability of finding three nucleons within the range of 
the correlator is sufficiently small to neglect three-body and higher-order contribu- 
tions to the cluster expansion. 

Different methods for the construction of optimal correlation functions were 
discussed in detail in Refs. [11,12]. Here we employ an energy minimisation in 
the two-body system using simple parametrisations for the central and tensor cor- 
relation functions. For the correlation functions R+(r) in the four possible (S, T) 
channels we will use one of the following forms: 



R A + (r) = r + a (r/fif exp[- exp(r/j8)] , 

R*(r) = r + a(l- exp[-(r/y)"]) exp[- exp(r/j0)] . 



The tensor correlations functions &(r) for the two 5 = 1 channels are parametrised 
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Table 1 

Parameters of the optimal central and tensor correlation functions for the AVI 8 potential as 

determined in Ref. [12]. 

Central correlation functions R+(r) 

(S,T) Param. a [fm] /? [fm] y [fm] 77 

(0,1) A 1.379 0.8854 — 0.3723 

(1.0) A 1.296 0.8488 0.4187 
(0,0) B 0.76554 1.272 0.4243 1 

(1.1) B 0.57947 1.3736 0.1868 1 







Tensor correlation functions &(r) 






(S,T) 


Param. 


a 


yS[fm] 


r [fm] 


V 


(l,0)a 


A 


530.38 


1.298 


1000 


1 


(1,0)7 


A 


0.383555 


2.665 


0.4879 


1 


(1,1) 


B 


-0.023686 


1.685 


0.8648 


1 



by 



#V) = a(l -exp[-(r/y)*])exp[-exp(r/jS)] , 
# B (r) = a(l -expKr/yfDexpt-r/jg] . 



(40) 



; 



The uncorrelated two-body wave functions are the free zero-energy solutions of the 
two-body Schrodinger equation with the lowest orbital angular momentum consis- 
tent with antisymmetry for given (S, T). The energy minimisation is performed for 
each (S, T) channel separately. The range of the tensor correlation function in the 
(S, T) = (1,0) channel is controlled by using the following integral constraint for 
the variation 

, „ 1 0.1 fm 3 ; a-correlator 

dr r 2 $(r) = \ ' . (41) 

I0.5fm ; y-correlator 

Depending on this constraint on the (5, T) = (1,0) tensor correlation functions, we 
will refer to the resulting sets of correlators as a- and y-correlators, respectively. 
The parameters of the correlation functions resulting from the energy minimisation 
for the AVI 8 potential are summarised in Table 1 . The details of their determination 
and the corresponding correlation functions for the Bonn A potential can be found 
in Ref. [12]. 

Plots of the optimal correlation functions for the AVI 8 potential are shown 
in Fig. 6. All central correlation functions R+(r) - r are of similar short range. 
Beyond r « 1 .5 fm for the even channels and r « 2 fm for the odd channels the 
shift R+(r) - r vanishes and the unitary correlation operator acts like the identity 
operator. In general, the maximum shift is smaller for the odd channels, since the 
uncorrelated L = 1 wave functions are already depleted at small r by the centrifugal 
barrier. The tensor correlation function §{r) in the (5, T) = (1,1) channel is very 
weak and does not lead to significant effects. For the dominant (S,T) = (1,0) 
channel, the optimal tensor correlators for the two different values of the constraint 
(41) are depicted. Only the a-correlator can be considered short-ranged. The y- 
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Fig. 6. Optimal central (a) and tensor correlation functions (b) for the AVI 8 potential. The 

different lines correspond to the different (S,T) channels: (0, 1) ( -), (1,0) ( ), 

(0,0) (- ),(1,1)( ). 

correlator has sizable contributions beyond 3 fm and one has to expect significant 
three-body contributions to the cluster expansion. 

An indirect estimate for higher-order contributions to the cluster expansion can 
be obtained by comparing the exact solution of a few-body problem using the cor- 
related Hamiltonian in two-body approximation with the solution based on the bare 
interaction. This has been done in the no-core shell model [8] for 4 He [18]. It turns 
out that for the ^-correlators the three- and four-body contributions to the correlated 
Hamiltonian are roughly 4 MeV. For the longer ranged y-correlators this contribu- 
tion grows significantly. 



4 Variational Ground State Calculations 

The correlated realistic interaction V UCO m provides a robust starting point for 
different approaches to tackle the many-body problem. The dominant short-range 
correlations are included explicitly though the unitary correlation operator so that 
simple low-momentum many-body spaces — not able to represent the short-range 
correlations themselves — suffice for a realistic description. 

4. 1 Variational Model — FMD States 

Here we will employ a variational model based on an extremely versatile para- 
metrisation for the many -body states which was developed in the framework of the 
Fermionic Molecular Dynamics (FMD) model [19,20,21,22]. The many -body trial 
state \Q) is given by a Slater determinant of single-particle states \q t ) 



\Q)=mqx)® |?2>®---® Wa))- 



(42) 



The single-particle states \q) are parametrised by Gaussian wave packets with vari- 
able spin orientation and fixed isospin. In general a superposition of several wave 
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Fig. 7. Energy expectation value for different nuclei calculated with the uncorrected Hamil- 
tonian (left), the central-correlated Hamiltonian (centre) and the fully correlated Hamilto- 
nian (right) in two-body approximation using the a-correlators. 

packets can be used to represent the single-nucleon state 



\q) = 2_jC v \a v , b v ) <8> \% v ) ® \m t ) 



(43) 



The Gaussian wave packet in configuration space is parametrised in terms of a 
complex width parameter a v and a complex vector b v = p v + ia v 7T v , where p v is the 
mean position and n v the mean momentum of the wave packet, 



(x\Oy, b v ) 



ex P[ " ~^^\ 



(44) 



In the simplest case each single particle- state contains 10 independent variational 
parameters. These are determined by a large-scale numerical minimisation of the 
expectation value of the correlated Hamiltonian 



E = 



(Q|H^-T cm |0 

(Q\Q) 



mm. 



(45) 



The operator T cm of the centre of mass kinetic energy is explicitly subtracted to 
eliminate the centre of mass contribution to the energy. 

The antisymmetrised Gaussian trial states are extremely flexible and allow for 
a multitude of intrinsic structures. Shell-model type states as well as states with 
strong ^-clustering can be described on the same footing. One should stress that in 
the following calculations none of these structures is put in by hand, but that they 
emerge naturally from the energy minimisation. A tremendous technical advantage 
is that all necessary matrix elements, e.g., for the various terms of the correlated 
interaction, can be calculated analytically. This makes large-scale variational cal- 
culations up to mass numbers A « 60 possible. Details on the matrix elements and 
the implementation are given in Ref. [23]. 

We use this variational model to illustrate the effect of the central and tensor 
correlations on the ground-state energy of different nuclei. To this end we per- 
form the energy minimisation with the fully correlated interaction (38) using the 
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a-correlators. With the resulting states we then compute the expectation values of 
the central-correlated Hamiltonian (without tensor correlations) and of the uncor- 
rected Hamiltonian. The results for 4 He, 16 0, 40 Ca, and 48 Ca are summarised in 
Fig. 7. The energy expectation value calculated with the uncorrelated Hamiltonian 
is large and positive, i.e., the system is not bound. As discussed in Sec. 1 the Slater 
determinant (42) is not capable of describing short-range central and tensor cor- 
relations. This entails that the repulsive core of the central interaction generates a 
large positive energy contribution and that the expectation value of the tensor part 
is practically zero. 

The proper inclusion of the short-range central correlations lowers the energy 
significantly, but still does not lead to bound nuclei since the attractive contribu- 
tion of the tensor interaction is missing. This is a striking demonstration of the 
importance of the tensor interaction and the associated correlations. Only after the 
inclusion of both, central and tensor correlation operators, we obtain bound ground- 
states. 

We point out that in general all tensor components of the correlated interaction 
Vucom detailed in (38) yield negligible expectation values with a single Slater de- 
terminant and can be omitted in the following. The tensor interaction enters only 
though the central, the 1 2 , and the (1-s) terms generated by the correlation procedure 
(cf. Sec. 3.2). For the sake of simplicity, the 1 2 (1 • s) contribution in (38) is replaced 
by 2(1 • s) and added to the conventional (1 • s) term. Thus two-body states with L = 
and L = 1 are treated correctly and contributions from L > 2 components, which 
are suppressed due to the centrifugal barrier, are approximated. 

4.2 Long-Range Correlations and Three-Body Interactions 

In the framework of the FMD variational model we calculate the ground state 
energies and the charge radii for a few nuclei using the correlated AVI 8 potential 
with the short-ranged a-correlators. Table 2 summarises the results and compares 
them to experimental data. The magnitude of the binding energies and the charge 
radii for larger nuclei are significantly smaller than the experimental values. There 
are two major reasons for these systematic deviations: long-range correlations and 
genuine three-body forces. 

Table 2 

Variational ground state energies and charge radii for the correlated AVI 8 potential with 

the a- and y-correlators in comparison with experiment. 





AYlSa 




AV18y 




Experiment 


Nucleus 


E/A [MeV] 


r C h IMI 


E/A [MeV] 


r ch [fm] 


El A [MeV] r ch [fm] 


4 He 


-4.18 


1.57 


-6.99 


1.51 


-7.07 1.68 


16q 


-4.07 


2.33 


-7.40 


2.25 


-7.98 2.71 


28 Si 


-3.45 


2.72 


-6.68 


2.66 


-8.45 3.12 


40 Ca 


-5.22 


2.87 


-8.19 


2.89 


-8.55 3.48 


48 Ca 


-5.14 


2.87 


-7.87 


2.93 


-8.67 3.47 
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Long-range correlations. Due to the restriction of the range of the tensor cor- 
relator in the (S,T) = (1,0) channel, long-range tensor correlations are not ad- 
equately described by the correlator. Remember, the range constraint is used in 
order to ensure the validity of the two-body approximation (cf. Sec. 3.5). Ideally, 
the missing long-range correlations would have to be accounted for by the model- 
space. However, in the case of the simple variational model, the single Slater de- 
terminant is not capable of describing long-range tensor correlations. In principle 
one could enlarge the model-space such that the long-range attraction provided by 
the tensor interaction can be exploited. This would lead to an increase in binding 
energy and radius but is presently not feasible, at least for larger nuclei. 

Alternatively one could resort to long-range tensor correlators at the price of 
including contributions of higher orders of the cluster expansion. A step into this 
direction is the use of the longer-ranged y-correlator in the (5, T) = (1,0) channel. 
The binding energies and charge radii resulting in two-body approximation are 
also given in Tab. 2. The binding energy increases significantly compared to the 
a-correlators. Part of this increase is an artifact of the two-body approximation: 
The three-body contributions, which are neglected here, would provide an effective 
repulsion and compensate most of the gain in binding energy. Hence, a treatment of 
long-range tensor correlations by the Unitary Correlations Operator necessitates a 
consistent inclusion of higher-order cluster contributions — a technically extremely 
challenging task. 

Genuine three-body forces. As quasi-exact Green's function Monte Carlo cal- 
culations for small systems show [5], realistic two-body potentials alone do not 
generate sufficient binding to reproduce experimental data. This can be remedied 
by introducing a phenomenological three-nucleon force which is adjusted to ground 
states and low-lying excitations [9]. Promising developments in effective field the- 
ories [10] might lead to realistic three-body forces which go beyond the rather 
phenomenological three-body terms considered so far in ab initio calculations. 



4.3 Phenomenological Correction 

The inclusion of three body terms, a genuine three-body force or three-body 
contributions of the cluster expansion, poses an enormous computational challenge. 
At the moment this is not feasible for the full range of particle numbers envisioned. 
We, therefore, resort to a pragmatic approach and employ a momentum-dependent 
two-body force to simulate the missing three-body terms and long-range correla- 
tions. The generic contribution of a momentum-dependent two-body interaction to 
the energy per particle in nuclear matter is similar to a local three-body force. 

The structure of the phenomenological correction should be as simple as possi- 
ble. For the single Slater determinant employed so far, only the local and momentum- 
dependent central terms and the spin-orbit term have sizable contributions. There- 
fore, a sensible ansatz for the correction is a sum of these three terms 

Av = Av c (r) + q Av w (r) q + Av^(r) (1 • s) . (46) 
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The symmetric form of the momentum-dependent part is chosen for convenience 
and can easily be transformed to the form appearing in (38) or (19). In order to keep 
the correction simple we consider only spin-isospin independent, i.e., Wigner-type, 
forces. Each of the radial dependencies Av(r) is given by a single Gauss function 

Av(r) = yexp[-r 2 /(2*)] (47) 

with strength y and range parameter k. The adjustment of the parameters of the 
three radial dependencies is done in two steps: First, the strengths and ranges of 
the central corrections are adjusted such that binding energies and charge radii of 
the doubly magic nuclei 4 He, 16 0, and 40 Ca are in agreement with experiment. 
Notice that there are only four parameters available to fit six quantities. Second, the 
parameters of the attractive spin-orbit correction are chosen such that the binding 
energies of 24 and 48 Ca are consistent with experiment. The resulting correction 
parameters for the correlated AVI 8 potential with the short-ranged a-correlators 
are 

Av c '(r) : y = -7.261 MeV, k = 2.75 fm 2 , 

Av w (r) : y = +14.05 MeV fm 2 , k = 2.5 fm 2 , (48) 

Av l \r) : y = -2.7 MeV, k = 3.0 fm 2 . 

It turns out that each of the three terms of the correction plays a different and 
physically quite intuitive role. The repulsive momentum-dependent term is respon- 
sible for the reduction of the central densities and the increase of the charge radii. 
The attractive local term provides the bulk of the missing binding energy. The at- 
tractive spin-orbit term gives additional binding especially for nuclei far off stabil- 
ity. Hence, (46) can be considered the simplest possible form of correction capable 
of generating all the necessary effects. 

One of the effects which is absorbed in the phenomenological correction is 
the lack of long-range tensor correlations. From this it is already clear that the 
correction depends on the available model space. For an extended model space, 
which is capable of describing part of the missing correlations by itself, the phe- 
nomenological correction will be different. The parameter set (48) is adapted to 
a model space of states consisting of a single Slater determinant as used in our 
variational approach or in Hartree-Fock calculations. If one allows for superposi- 
tions of Slater determinants, e.g., in the framework of variation after projection or 
multi-configuration calculations, the correction has to be readjusted. 

4.4 Binding Energies and Charge Radii 

Using the a-correlated AVI 8 potential with the phenomenological correction 
(46) we perform variational ground state calculations in the mass region A < 60. 
Figure 8 depicts the part of the nuclear chart accessible to the variational model 
based on the FMD trial state (42). The shadings indicate the deviation of the energy 
expectation value (45) from the experimental binding energy per nucleon [24]. The 
main chart was obtained using a single Gaussian wave packet to describe the single- 
particle states. The agreement with experiment of the absolute binding energies 
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Fig. 8. Chart of nuclei accessible to the variational model. The shading shows the differ- 
ence between the variational and the experimental binding energy. The main chart (a) was 
obtained using a single Gaussian wave packet per nucleon. The inset (b) shows results with 
two Gaussians per nucleon. 

around the magic numbers and for larger nuclei is quite good (bright shades). The 
largest deviations appear for p- and sd-shell nuclei away from the shell closures. 
One reason for these deviations is the insufficient flexibility of the single-particle 
trial states to simulate, e.g., the long-range behaviour of the states especially for 
light nuclei. Enhancing the flexibility by using a superposition of two Gaussians 
to parametrise the single-particle states leads to a significant improvement, as the 
small chart in Fig. 8 demonstrates. 

A more quantitative view is given in Fig. 9, where the binding energy devia- 
tions (E - E exp )/A and the charge radii r c h are shown for selected stable isotopes. 
The calculated charge radii are in excellent agreement with the experimental data 
[25] throughout the whole mass range. The binding energies are in nice agreement 
with experiment in the vicinity of the magic nuclei but start to deviate in between 
the shell closures, especially for p-shell nuclei. As already mentioned, part of this 
deviation is due to an insufficient flexibility of the single-particle trial states. A re- 
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Fig. 9. Difference between variational and experimental binding energy per nucleon (up- 
per panel) and charge radius (lower panel) for selected stable isotopes. All calculations 
are performed using the a-correlated AVI 8 potential plus phenomenological correction 

(46). Shown are results obtained with one ( ) and with two Gaussian wave-packets 

( ) for each single-particle trial state. 
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Fig. 10. Binding energy difference and charge radius for Calcium isotopes (cf. Fig. 9). De- 
picted are results based on the cK-correlated AVI 8 potential with the full phenomenological 
correction ( ) and with the central terms of the correction only ( ). 

finement of the trial states by using a superposition of two Gaussian wave-packets 
(dashed line) instead of one (solid line) leads to a noticeable improvement of the 
energies. The remaining discrepancy can be reduced further by angular momentum 
projection as discussed in Sec. 5. 

For larger nuclei, the trial state with one Gaussian per nucleon provides a good 
description. The binding energy differences and charge radii for the Calcium iso- 
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topes are depicted in Fig. 10. The binding energies are in very good agreement with 
experiment, the charge radii are slightly overestimated for the heavier isotopes but 
are generally in agreement with experiment. Figure 10 also illustrates the effect of 
the spin-orbit contribution in the phenomenological correction. The gray symbols 
show results of a calculation using the or-correlated AVI 8 potential with central 
correction but without the spin-orbit term. The spin-orbit term provides additional 
attraction for the N # Z isotopes and reduces the charge radii at the same time. The 
closed- shell N = Z nuclei are not affected. 

4.5 Intrinsic Single-Particle Density Distributions 

Since the variational calculation provides us with the full many-body wave 
function we can easily compute other physical quantities, e.g., the intrinsic single- 
particle density distribution 



P(x) = 2 (21 T^ m ,(x)Y mi , m ,(x) |0 . (49) 



m s ,m t 



Unlike for the two-particle density distribution, the short-range central and tensor 
correlations have only marginal influence on the diagonal elements of the one- 
particle density matrix [11,14]. Therefore, we restrict ourselves to the uncorrelated 
one-body densities for simplicity. 

Figure 1 1 shows three-dimensional illustrations of the intrinsic density distribu- 
tions for selected nuclei. Depicted is the iso-density surface forp(jt) = po/2, where 
p = 0.17 fm is the nuclear matter saturation density. In addition one octant is 
removed to visualise the interior density distribution. For the doubly magic nuclei 
16 0, 40 Ca, and 48 Ca the variational calculation leads to spherical symmetric density 
distributions in accord with shell-model-type states. For the p- and sd-shell nuclei 
7 Li, 9 Be, 12 C, 20 Ne, and 26 Mg the one-body density distributions exhibit pronounced 
intrinsic deformations. The density distribution of 9 Be reveals a clear two-a sub- 
structure; 20 Ne shows a toroidal belt similar to a 12 C nucleus supplemented with two 
a-clusters forming end-caps. For 26 Mg only remnants of a-clustering are visible in 
the rather smooth density distribution with a multipolar deformation. 

These plots showcase the flexibility of the FMD states. A Slater determinant of 
Gaussian single-particle states is capable of describing shell-model wave functions 
as well as states with strong a-clustering. We should stress that the structures shown 
in Fig. 1 1 are generated by a subtle interplay between the different terms of the 
realistic potential. They are not imposed through constraining the trial state like in 
a-cluster models. 

In many cases, strong intrinsic deformations and a-clusters present in a N « Z 
nucleus dissolve gradually when neutrons are added (or removed). An example are 
the Neon isotopes depicted in Fig. 12. Starting from 20 Ne with its characteristic 
intrinsic deformation, the addition of neutrons washes out the a-clusters and even- 
tually leads to an almost spherical, shell-model-type density distribution for 26 Ne. 

Finally, in Fig. 13 we compare the radial charge density distributions for 12 C, 
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Fig. 1 1. Intrinsic total density distributions of several stable nuclei. Shown is the iso-density 
surface for p(x) = po/2 with one octant cut away to reveal the density distribution in the 
centre. The colour coding on the section gives the density according to the colour bar at the 
right. The background grid has a mesh size of 1 fm x 1 fm. 




Fig. 12. Intrinsic single-particle density distribution p(x) for different Neon isotopes (cf. 
Fig. 11). 

16 0, 40 Ca, and 48 Ca with a model-independent analysis of electron scattering data 
[25]. For this comparison, the exact centre of mass correction and the proton and 
neutron form factors are included in the theoretical curves. The calculated density 
profiles are in good agreement with the experimental analysis. The surface struc- 
ture is reproduced extremely well. In the interior the modulations of the calculated 
density distribution are slightly more pronounced. This behaviour is common to 
most models based on a single Slater determinant for the ground state; one possi- 
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Fig. 13. Radial charge distribution p c ^{x) for different nuclei obtained from the varia- 
tional wave function ( ) in comparison with the experimental charge distribution 

( ) resulting from a model-independent analysis [25]. 

bility to improve upon this are multi-configuration calculations. Furthermore, the 
uncertainties in the experimental determination of the density profiles are largest in 
this region. 



5 Angular Momentum Projection 

5. 1 Projection after Variation 

Many of the density distributions presented in Sec. 4.5 reveal pronounced in- 
trinsic deformations. Evidently the states obtained by energy minimisation for the 
trial state (42) are not necessarily angular momentum eigenstates as it is the case 
for the eigenstates of the Hamiltonian — they have to be interpreted as intrinsic 
states. One can extract angular momentum eigenstates I^Vjm) from the intrinsic 
state 12) by a standard angular momentum projection technique [26]. The angular 
momentum projected states are defined through 



\^jm) = Y,S J kKk\Q) 



(50) 



The generalised projection operator P J MK is given by 

27+1 



MK 



Sn 2 



fdcoD 



J * K (a>) R(a>) 



(51) 
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where D J MK (a>) are the Wigner D-functions and R(o») is the unitary rotation operator 
expressed in terms of the three Euler angles a> = (a,fi, y): 

R(aS) = exp(-iaJ z ) exp(-i/?J v ) exp(-iyJJ . (52) 

The coefficients g J K are formally determined by minimising the energy expectation 
value in the projected state I^jm) 

F _ C^jmI H - T cm YVjm) _ Yakk 1 8k Sk'^kk' ,,,> 

/_ ~ wjmWjm) ~~i:KK'8 J K *8> KK r 

where h J KK , = (Q\ (H - T cm )P^, \Q) and n J KK , = (Q\ P J KK , \Q). This leads to a gener- 
alised eigenvalue problem for the energies Ej of the projected states and the coef- 
ficients g J K : 

YjhKK'Si^EjYj^KK'SK'- (54) 

K' K' 

In this paper we employ the angular momentum projection in the framework 
of a projection after variation (PAV) calculation, i.e., we use the intrinsic state \Q) 
obtained by minimising the energy expectation value (45) and subsequently project 
it onto angular momentum eigenstates (50). 

This scheme is different from a variation after projection (VAP) calculation 
which is based on a minimisation of the energy expectation value (53) calculated 
with the projected states. The presence of strong intrinsic deformations is associ- 
ated with a kinetic energy contribution that makes such states unfavourable when 
minimising the intrinsic energy. This energy contribution, however, is removed by 
the projection procedure. Therefore, clustering and intrinsic deformation are gen- 
erally more pronounced in a VAP framework, where the projected energy is min- 
imised. In terms of the available model-space, a VAP calculation goes beyond the 
variation with a single Slater determinant. The VAP trial state is a specific super- 
position of rotated Slater determinants. As discussed in Sec. 4.3 this, in general, 
necessitates a readjustment of the phenomenological correction. First results of 
a variation after parity projection calculation are available. A full variation after 
angular momentum projection is computationally extremely involved and will be 
discussed elsewhere. 

The effect of the angular momentum projection on the ground state energy is 
illustrated in Fig. 14, where the intrinsic energies are compared to the energy of 
the lowest angular momentum eigenstate obtained in the PAV framework. The de- 
viation from the experimental binding energies is reduced by a factor 0.5 for the 
intrinsically deformed p- and sd-shell nuclei. The energy of spherical nuclei is not 
affected. The residual difference has to be accounted for by further extending the 
model space, e.g., in the framework of a multi-configuration calculation. 

5.2 Rotational Spectra 

Besides the ground state the PAV procedure provides us with a whole rotational 
spectrum based on the intrinsically deformed variational state. Figure 15 shows 
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Fig. 14. Energy deviation (E - E sxp )/A for the intrinsic variational states ( ) and for 

the angular momentum projected intrinsic states (projection after variation) ( ). For 

A < 28 single-particle states composed of two Gaussians are used, heavier isotopes are 
calculated with one Gaussian per nucleon. 

spectra obtained by angular momentum projection of the intrinsic states for 7 Li, 
9 Be, 12 C, and 20 Ne. To give an impression of the influence of the phenomenological 
correction to the a-correlated AVI 8 potential, three different calculations are shown 
for each nucleus. First, the a-correlated AVI 8 potential is used without correction 
term (left column in each panel). As discussed in Sec. 4.2 binding energies and 
charge radii are generally too small in this case. Nevertheless, the relative spectrum 
agrees very well with the experimental one. This is a remarkable result, considering 
that for the rotational band the moment of inertia and thus the radius of the intrinsic 
density distribution is crucial. 

Second, we perform the variation and projection with the a-correlated AVI 8 
potential including only the central part of the phenomenological correction (46). 
Compared to the significant influence of the correction on the binding energy and 
the radius of the ground state the effect on the spectrum is minor (cf . second column 
in each of the panels in Fig. 15). Generally the spectra are slightly compressed due 
to the increase of the radii induced by the central part of the phenomenological cor- 
rection. The inclusion of the spin-orbit part of the correction stretches the rotational 
bands again (see third columns in Fig. 15). 

Overall, the agreement with experimental data is very encouraging. The ground 
state rotational bands of 12 C and 20 Ne are reproduced very well. Further improve- 
ments are possible in the framework of multi-configuration calculations, where 
some of the missing states corresponding to collective modes can also be described 
[27]. 



6 Summary and Outlook 

We have combined two powerful tools for the treatment of the nuclear many- 
body problem: the Unitary Correlation Operator Method (UCOM) and the Fermionic 
Molecular Dynamics (FMD) model. 

The Unitary Correlation Operator Method (UCOM) provides a systematic way 
to derive, from realistic NN-potentials, effective interactions suitable for simple 
(low-momentum) model- spaces. Dominant central and tensor correlations induced 
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Fig. 15. Energy spectra resulting from the angular momentum projection of the variational 
ground states for 7 Li, 9 Be, 12 C, and 20 Ne. Shown are the excitation energies obtained using 
the a-correlated AVI 8 potential without correction (first column), with central correction 
only (second column), with central and spin-orbit collection (third column) in comparison 
to the experimental spectra (fourth column). The numbers given beneath the lowest level 
are the absolute ground state energies in [MeV]. 

by the NN-potential are described explicitly by an unitary transformation. By ap- 
plying the unitary correlation operator to a simple many-body state, e.g., a Slater 
determinant, the short-range correlations are imprinted into the state. Alternatively 
the Hamiltonian (and all other operators) can be transformed. The resulting corre- 
lated interaction Vucom possesses a closed operator representation, is phase-shift 
equivalent to the original potential, and leads to momentum- space matrix elements 
in accord with the Vi ow _£ matrix elements [12]. It provides a robust starting point for 
many -body methods relying on model spaces not capable of describing short-range 
correlations themselves, e.g., variational models with simple trial states, Hartree- 
Fock calculations, and multi-configuration shell-model. 

We have employed the correlated interaction Vucom m the variational frame- 
work of the Fermionic Molecular Dynamics (FMD) model. The many-body trial 
state is a Slater determinant of Gaussian single-particle states which prove to be 
extremely flexible. On the same footing, it allows for the description of spheri- 
cal shell-model-type wave functions as well as states with intrinsic deformations 
and a-clustering. At the same time, the approach is computationally quite efficient 
so that a treatment of nuclei up to mass number A < 60 is possible for realistic 
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two-body interactions with complex operator structure. The inclusion of three-body 
forces would substantially increase the computational cost and thus reduce the ac- 
cessible region of the nuclear chart. At the present stage, three-body interactions 
and long-range correlations are, therefore, simulated by a phenomenological two- 
body correction. Charge radii and charge distributions are in nice agreement with 
experiment. Around the magic numbers the variational energies agree well with 
the experimental binding energies. Away from the shell closures the intrinsic den- 
sity distributions exhibit strong deformations and a-clustering, which necessitates 
projection of the intrinsic states onto angular momentum eigenstates. Within a pro- 
jection after variation (PAV) framework we have achieved very promising results 
for ground state energies and rotational spectra. 

Natural next steps are variation after projection (VAP) and multi-configuration 
calculations. First results along these lines have been obtained using a Generator 
Coordinate Method to implement an approximate VAP scheme. Constraints on the 
multipole moments are used to generate different intrinsic states which serve as 
input for variation after projection (w.r.t. the generator coordinates, i.e., the mul- 
tipole moments) and multi-configuration calculations. VAP calculations for 16 
show, e.g., that a tetrahedron of a-clusters is energetically more favourable than 
the spherical shell-model-type state resulting from the unprojected variation [27]. 
Detailed results of these investigations will be discussed in a forthcoming publica- 
tion. 

Other future applications of the correlated realistic interaction Vucom are Hartree- 
Fock and RPA calculations. This provides a consistent scheme to perform nuclear 
structure calculations based on realistic NN-potentials also for large nuclei. 
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